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Effect of initial conditions on aging properties of the spin-glass state is studied for a single crystal 
Cu:Mn 1.5 at %. It is shown that memory of the initial state, created by the cooling process, 
remains strong on all experimental time scales. Scaling properties of two relaxation functions, the 
TRM and the IRM (with t w i — and t W 2 = t w ), are compared in detail. The TRM decay exhibits 
the well-known subaging behavior, with ^ < 1 and > t w . The IRM relaxation demonstrates the 
superaging behavior, with fi > 1 and < t w . It is shown that an average over different initial 
conditions leads to systematic improvement in t/t w scaling. An effective barrier, describing influence 
of the initial state on spin-glass relaxation, is found to be almost independent of temperature. These 
results suggest that departures from full t/t m scaling, observed in aging experiments, are largely due 
to cooling effects. 
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I. INTRODUCTION 

Aging phenomena in rSpin glasses have been studied 
for almost twenty years.EI Despite considerable progress, 
some problems remain open. One of them is a problem of 
t/t w scaling of time dependent quantitiesH This issue is 
a touchstone of our understanding of spin-glass dynam- 
ics. Only when all departures from full t/t w scaling are 
accounted for, both theoretically and experimentally, can 
one say that the aging phenomena are well understood. 

First experimental results on t/t w scaling of the ther- 
moremanent magnetization (TRM) were obtained soon 
after the aging effects had been discovered. They demon- 
strated that there are small but systematic deviations 
from full t/t w scaling in the aging regime.ErQ It was shown 
that the TRM curves for different waiting times can be 
successfully superimposed if a power of the waiting time, 
i^,, with ji < 1, is used in the analysis instead of the 
actual t w . This phenomenon is often referred to as "sub- 
aging". The apparent age of the spin-glass state, de- 
termined from the scaling of the TRM curves, increases 
slower than the waiting time t w . It was also showntffl that 
the quasiequilibrium decay at short observation times 
should be taken into account properly in order to achieve 
good scaling over a wide time range. However, full t/t w 
scaling in the aging regime has never been derived from 
experimental data. Deviations of /i from unity are small, 
but persistent. The physical meaning of this scaling pa- 
rameter remains unclear eyeti though different interpreta- 
tions have been proposed.Em All this suggests that there 
may be some additional reasons for the lack of scaling, 
which have not been taken into account yet. 

On the theoretical side, the situation is mor* ambigu- 
ous. Phcnomcnological phase-space models,!™ describ- 
ing aging as thermally activated hopping over free-energy 
barriers, suggest full t/t w scaling in the absence of finite- 



size effects. In the mean-field dynamics, however, two 
different cases are distinguished. El Those mean-field mod- 
els, where one level of replica-symmetry breaking is exact 
(like the spherical p-spin model), have only one time scale 
in the aging regime, and t/t w scaling is predicted.lla Mod- 
els with continuous replica-symmetry breaking (like the 
Sherrington-Kirkpatrick modeU_ are characterized by an 
infinite number of time scales .Eil If decay of the correla- 
tion function in the aging regime is viewed as a sequence 
of infinitesimal steps, ithen each step takes much longer 
than the previous one.E3 Therefore, no t-Jj-^ scaling is ex- 
pected in these models. It is suggestecoid that full t/t w 
scaling for all times in the aging regime would rule out the 
continuous replica-symmetry breaking scenario. Thus, a 
reliable experimental evidence of presence (or absence) 
of t/t w scaling may help determine validity of different 
theoretical pictures. 

Numerical studies of spin-glass dynamics reveal an in- 
teresting fact. The short-range 3-dimensional Edwards* 
Anderson (EA) model exhibits good t/t w scaling,E£tEfl 
while there is no such scaling in .the .infinite-range 
Sherrington-Kirkpatrick (SK) modcl.E3£jO However, va- 
lidity of the EA model for description of rcaJ spin-glasses 
has been questioned. Recent simulationsliia have shown 
no trace of rejuvenation effects, observed in temperature- 
change experiments.113 It was argued that either time 
scales, reached in simulations, are too short, or the 
Edwards- Anderson model is not a good model for real 
spin glasses.Ea Further studies of both t/t w scaling and 
temperature-variation effects are needed to find out 
which model is better suited for theoretical description 
of spin-glass dynamics. 

The present paper is devoted to a study of departures 
from full t/t w scaling in spin-glass aging experiments. We 
hypothesize that one of possible reasons for the lack of 
such scaling is the influence of the cooling process. Our 
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line of argument is the following. A typical spin-glass re- 
laxation experiment includes cooling from above the glass 
temperatureEa Approach of the measurement tempera- 
ture is necessarily slow with possible oscillations, because 
temperature. Has. Xq be stabilized. Temperature-cycling 
experimentsjHduj'E^I as, well as aging experiments with dif- 
ferent cooling rates,Eil have shown that thermal history 
near the measurement temperature has a profound effect 
on the subsequent spin-glass behavior. Any aging exper- 
iment is, therefore, a temperature-variation experiment, 
and measured relaxation properties are determined by 
both cooling and waiting time effects. As the waiting 
time increases, the influence of the initial condition, set 
by the cooling process, diminishes. This may lead to a 
systematic deviation from full t/t w scaling, compatible 
with the experimentally observable behavior. In order 
to extract relaxation curves that exhibit good t/t w scal- 
ing, one has to average experimental results over different 
initial conditions. 

The paper is organized as follows. The next Section 
discusses general properties of spin-glass relaxation and 
methods for analysis of t/t w scaling. In Sec. III. A, fea- 
tures of the cooling process and properties of the TRM 
decay for t w = are studied. Sec. III.B discusses two dif- 
ferent types of deviations from t/t w scaling. In Sec. III.C, 
properties of the TRM and the IRM are analyzed in de- 
tail. The effect of initial conditions on scaling properties 
of measured relaxation is discussed in Sec. III.D. The last 
Section summarizes our arguments. 



II. THEORETICAL BACKGROUND 
A. Linear response in spin glasses 

In a typical TRM experiment, a spin-glass sample is 
cooled down from above the glass temperature T g to a 
measurement temperature T < T g in the presence of a 
small magnetic field h. It is then kept at the measure- 
ment temperature for some time t w , which is called wait- 
ing time. After that, the field h is cut to zero, and a 
decay of the thermoremanent magnetization (TRM) is 
measured as a function of observation time t, elapsed af- 
ter the field change. This decay depends on the waiting 
time t w — a phenomenon called aging. The total age of 
the system is t + t w . In a more complex IRM proto- 
col, the sample is cooled down at zero magnetic field and 
kept for a waiting time t w \. Then a small magnetic field 
is turned on, and, after an additional waiting time t w2 , is 
turned off again. A subsequent decay of the isothermal 
remanent magnetization (IRM) is observed. It depends 
on both waiting times. 

In order to describe time evolution of the system, the 
autocorrelation function, C(ti,t 2 ), is introduced: 

JV 

C(h,t 2 ) = (l/N)J2(S t (h)S l (t 2 )) . (1) 
i=l 



It depends on two times, t\ and t 2 , measured from the 
end of the cooling process. Response of the system at 
time t\ to an instantaneous field h, present at time t 2 , is 
described by the response function: 

N 

R(t 1 ,t 2 ) = (l/N)Y,5(S i (t 1 ))/5h(t 2 ) . (2) 

i=l 

In thermal equilibrium, both functions depend only on 
the time difference, t\ — t 2 , and they are related by the 
fluctuation-dissipation theorem. If the system is out of 
equilibrium, a generalization of this theorem is expected 
to holdfl 

R(t 1 ,t 2 )=f3X[C]dC(t 1) t 2 )/dt 2 . (3) 

Here /3 = l/fc^T, and X[C] is the fluctuation-dissipation 
ratio, which is equal to unity in equilibrium. It is 
suggestedQ that, for long waiting times, X depends on 
its time arguments only through the correlation function, 
i.e. X(t 1 ,t 2 ) = X[C(t 1 ,t 2 )]. 

Magnetic susceptibility, measured in spin-glass exper- 
iments, is the integrated response. If the field is present 
during the waiting time, t w , the susceptibility, measured 
at the observation time, t, is given by the following ex- 
pression: 

X(t + t w ,t w )= f R(t + t w ,t')dt' . (4) 
Jo 

Using Eq. (3) for the response function and introducing a 
function Y[C] through a relation (3X[C] = dY/dC$ one 
obtains the formula: 

x(t + t w ,t w )=Y[C(t + t w ,t w )]-Y[C(t + t w ,0)] . (5) 

The first term on the right is. identified with the in- 
dependent TRM susceptibility.cl The second term corre- 
sponds to the TRM decay after zero waiting time. To 
avoid confusion, we shall refer to it as ZTRM, and drop 
the second argument. The function on the left is the 
IRM susceptibility for t w \ = and t w2 = t w . In what 
follows, the IRM will always denote the isothermal rema- 
nent magnetization with t w \ = 0. It is characterized by 
the same waiting time, t w , as the TRM. We shall also use 
a corresponding observation time instead of a total age as 
the first argument of each function. For the ZTRM, the 
total age is equal to the observation time. Thus, Eq. (5) 
can be rewritten in the following form: 

— XZTR,Ai(t + tw) ■ (6) 

The last formula expresses the well-known principle of 
superposition, which has_been verified experimentally in 
the case of spin glasses. cB 

It follows from Eq. (6) that the TRM relaxation is a 
superposition of two decays. The IRM is a response, as- 
sociated with the waiting time. The ZTRM is a response, 
related to the cooling process. However, the ZTRM it- 
self cannot be treated as a linear response because of the 
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changing temperature. This means that the thermal his- 
tory cannot be taken into account by simple addition of 
the cooling time to the waiting time. 

If memory of the initial state is strong, there will be 
no one-to-one correspondence between the i^-dependent 
correlation and i^-dependent linear response. The TRM 
depends on the correlation C(t + t w ,t w ), but it includes 
response of the initial state. The IRM is a linear re- 
sponse, associated with the waiting time only, but it de- 
pends on the correlation with the initial state. Therefore, 
neither the TRM, nor the IRM can be called the "true" 
independent response. 

The mean-field theory of aging phenomena relies on 
the assumption of weak long-term memory. According 
to this assumption, the system responds to its past in an 
averaged way, and its memory of any finite time interval 
in the past is weak .EH This implies that the lower integra- 
tion limit, t' = 0, in Eq. (4) can be replaced by any finite 
t' = to, as long as both t and t w go to infinity.EJ Thus, 
there should be no difference between the TRM and the 
IRM in the long-time limit. Within this asymptotic, s 
proach, the IRM is sometimes referred to as TRM.E3E 

The linear response theory predicts that the measured 
susceptibility, Eq. (4), depends on the waiting time t w . 
However, it gives no predictions regarding t/t w scaling. 
Eq. (6) suggests that the TRM and the IRM, though 
characterized by the same waiting time, cannot have the 
same scaling properties. This is because they differ by 
the ZTRM, which is a one-time quantity without a char- 
acteristic time scale. Two exceptional cases are possible. 
First, the ZTRM decays so rapidly, that the last term in 
Eq. (6) can be neglected for sufficiently large Lp.. This 
case is usually considered in theoretical studiesjja which 
assume that memory of the initial condition is lost af- 
ter very long waiting times, i.e. C(t + t w ,0) — > as 
t w — > oo.Ej Second, the ZTRM changes so slowly, that 
XzTRKi(t + tw) for long enough t w can be treated as a 
non-zero constant. This approach is similar to the treat- 
ment of the field-cooled susceptibility, which is also a 
one-time quantity. We show in Sec. III.B that neither of 
these two conditions holds in spin-glass experiments. 



B. Analysis of t/t w scaling 

Spin-glass dynamics is characterized by at least two 
time scales. The microscopic attempt time, To, is associ- 
ated with the quasiequilibrium decay at short observation 
times. The waiting time, t w , determines properties of the 
nonequilibrium relaxation at long times. Both types of 
spin-glass behavior have to be taken into account when 
t/t w scaling of relaxation curves is analyzed. 

Numerical studies of the off-equilibrium dynamics in 
the 3D Edwards- Anderson modefij suggest that the au- 
tocorrelation function, C(t,t w ), can be well approxi- 
mated by a product of two functions: 



The waiting time independent factor t~ a represents the 
slow (on the logarithmic scale) quasiequilibrium decay at 
t -C t w . The function $(t/t w ) (which is approximately 
constant for t -C t w ) describes the faster nonequilibrium 
relaxation at longer observation times. Of course, both 
a and <j>(x) depend on temperature. 

A similar multiplicative ansatzE3 has been successfully 
used for scaling experimental relaxation curvesnu 



XTR M {t,t w ) cxt- a F(t e / 



(8) 



Two features distinguish Eq. (8) from Eq. (7). First, 
an effective time t e (usually denoted by A) is introduced. 
The age of the system increases with the observation time 
as t w + 1, and the time/age ratio decreases. In order to 
allow description of the relaxation by the same age t w , the 
effective time should increase slower than the observation 
time t. This means that dt e /t w = dt/(t + t w ) and t e = 
t w hi(l+t/t w ). Second, possible deviations from full t/t w 
scaling in the aging regime are taken into account by the 
parameter /i. In this case, dt e /t£ = dt/(t + t w )^, and the 
effective time is 



t e = t w [{l + t/t w ) 



i]/(i-m) 



(9) 



At short observation times, t <C t w , the effective time is 
equivalent to the observation time: t e w t. 

The /^-scaling approach is very useful for studying de- 
partures from t/t w scaling, no matter whether /j, itself 
has a clear physical meaning or not. We use this method 
in the present paper and determine \i from juxtaposition 
of relaxation curves, plotted versus t e /t£. 

A different method for separating the quasiequilib- 
rium and— aging regimes has also been successfully 
cmployed.EI I^|is, ^spired by dynamical solution of mean- 
field models 
is considered 




The following additive representation 



XTR M {t,t w ) OC A{t/T Q y a + f(t e /t%) 



(10) 



C(t,t w ) cx r a $(t/t v 



(7) 



This approach yields better t/t w scaling £ji_. closer to 
unity) than the previously discussed methodE'cl It should 
be noted, however, that derivation of the last formulacl 
employs an assumption that C(t + t w ,0) = 0. Therefore, 
Eq. (10) is exact only asymptotically (large t w ), when 
long-term memory is weak. Nevertheless, this method is 
justified, because it gives good results on experimental 
time scales. 



III. EXPERIMENTAL RESULTS AND 
ANALYSIS 

The purpose of this paper is to study how different 
initial conditions affect aging phenomena in spin-glass 
experiments. By the initial condition we mean history of 
the spin-glass state before the waiting time begins. We 
do not attempt to prove existence or absence of full t/t w 
scaling. We try to understand, on the qualitative level, 
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where observable departures from the full scaling in the 
aging regime come from. 

All experiments were performed on a single crystal of 
Cu:Mn 1.5 at %, a typical Heisenberg spin glass with a 
glass temperature of about 15.2 K. The single crystal 
was used to avoid possible complications due to finite- 
size effects. Another advantage of this sample is its 
high de Almeida-Thouless (AT) critical line.E3 For exam- 
ple, the AT field at T/T g = 0.87, the highest measure- 
ment temperature in our experiments, is about 600 Oe. 
This enables us to use a relatively large field change, 
AH = 10 Oe, and still work well within the linear re- 
sponse regime. A commercial Quantum Design SQUID 
magnetometer was used for all the measurements. This 
equipment has been optimized for precision and repro- 
ducibility, rather than speed. Thus, cooling is relatively 
slow, and this allows us to study cooling effects in detail. 



A. Cooling process and ZTRM 

Spin-glass behavior is characterized by memory effects: 
properties of measured spin-glass relaxation depend on 
history of a sample below T g . The cooling process is an 
integral part of this history. It has been shown that, in 
real spin glasses, details of the experimental protocol well 
above the measurement temperature, T, do not affect the 
measured relaxationHJ However, the thermal history in 
the immediate vicinity of the measurement temperature 
(say, ST < 0.5 K) has a strong impact on the observed 
spin- glass behavior. This fact suggests that the cooling 
and waiting time effects cannot be separated and should 
be studied simultaneously. 

A typical cooling protocol is exhibited in the inset of 
Fig. 1. The temperature drops rapidly to well below 
the measurement temperature, than rises to T + ST, and 
slowly approaches T from above. By definition, the ex- 
perimental time starts when the measurement tempera- 
ture, T, is finally reached. 

The approach of the measurement temperature is 
shown in the main body of Fig. 1. For relatively high 
temperatures, ST « 0.1 K. For the lowest tempera- 
ture, however, ST ~ 0.3 K. Comparison with results 
of temperature-cycling experiments^!] suggests that the 
initial undercool is not very important, because of its 
large (several Kelvin) magnitude. The subsequent over- 
shoot, however, may be expected to play a significant role 
in spin-glass dynamics, and cannot be neglected in our 
experiments. It should be noted that Fig. 1 exhibits the 
temperature curves for helium gas, used as a heat transfer 
agent within a sample chamber. The actual temperature 
of the spin-glass sample lags behind. This means that, in 
reality, the undercool is probably less pronounced, and 
the sample spends more time above the measurement 
temperature. 

It is instructive to compare the cooling procedure in 
Fig. 1 with an experimental protocol including a pos- 
itive temperature cycle.EJ In this protocol, the sample 
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FIG. 1: Approach of the measurement temperature, T, dur- 
ing the cooling process. The temperature curves T'(t), from 
top to bottom, are for T/T g = 0.3, 0.4, and 0.77, respectively. 
The inset shows the whole cooling process for T/T g = 0.77. 



is field-cooled to the measurement temperature, T, and 
kept for a long waiting time t w \. Then it is heated up to 
T + ST and annealed for a short time t w %. After that, it 
is cooled down and kept at the initial temperature T for 
additional time t w3 . Then the field is cut to zero, and the 
TRM relaxation is measured. It turns out that this re- 
laxation follows a reference curve with t w = t W 3 at short 
observation times, then breaks away and moves towards 
a reference curve with t w — t w ±. Therefore, the spin- 
glass state initially exhibits memory of only those events 
that happened after the temperature cycle. But, as time 
progresses, it begins to recall its earlier history. Another 
relevant analogy is an experimental protocol with a neg- 
ative temperature shift after initial waiting.Ej In both 
experiments, the relaxation curve after the temperature 
change cannot be merged with any of the reference curves 
measured at a constant temperature. 

The main idea of this paper is that any spin-glass aging 
experiment is, at the same time, a temperature-cycling 
experiment. All effects, which manifest themselves in 
temperature-cycling experiments, also appear in aging 
experiments. Our reasoning is based on thcphi^archical 
phase-space picture of spin-glass dynamics Let us 
consider a typical cooling process and imagine that the 
system spends some short time at temperature T + ST. 
Several metastable states, separated by free-energy bar- 
riers, become populated during that time. As tempera- 
ture is lowered, these states split in a hierarchical fash- 
ion and produce new metastable states with new barri- 
ers. The old barriers grow steeply, and, for any value 
of temperature, there exist barriers, diverging at this 
temperatureBj Therefore, when the measurement tem- 
perature T is finally reached, there is a large number of 
metastable states, separated by barriers of all heights. In 
other words, the cooling process creates an initial state 
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with a broad spectrum of relaxation times. This situ- 
ation corresponds to "existence of domains of all sizes 
within the initial condition" in the real-space picture.ta 

If a waiting time, t w , follows the cooling, properties 
of the subsequent relaxation will be similar to those in 
the positive-temperature-cycle experiment. At short ob- 
servation times, the relaxation will be governed by t w . 
At longer observation times, it will break away and slow 
down, because the long-time metastable states with high 
barriers, created by the cooling process, will come into 
play. This effect will be more pronounced after shorter 
waiting times. 

It is important to mention that the results of the 
temperature-variation experiments can be explained 
within the phcnomenological real-space picture, if 
a strong .separation of time and.— length scales is 
postulatedBj The experimental data,E3 that suggest di- 
vergence of free-energy barriers at any temperature be- 
low T g , can be reanalyzed in terms of barriers, which, 
vanish at T g , and remain finite at lower temperatures.ca 
The barriers, nevertheless, grow with decreasing temper- 
ature. This fact is sufficient for understanding results of 
the present paper. 

Properties of the initial state, produced by the cooling 
process, can be studied by measuring a TRM decay for 
t w = 0, which we call ZTRM. The sample is cooled down 
in the presence of a small field, H — 10 Oe. When the 
temperature is stabilized, the field is cut to zero, and 
decay of the ZTRM is recorded. 

Fig. 2 shows the logarithmic relaxation rates of this de- 
cay, corresponding to different temperatures. Each curve 
has a peak, and we call a position of this peak the "ef- 
fective cooling time" , t e p . Each magnetization curve is 
normalized by one at t = ff> before differentiation is 
performed. It can be seen from Fig. 2 that the peaks 
become much broader as temperature decreases. This 
means that the spectrum of relaxation times broadens as 
well, and it is not dominated by a single time scale. 

Even though the effective cooling time, tf", corre- 
sponds to the maximum of the relaxation rate, it cannot 
be treated as a regular waiting time t w . We have tried 
to scale each ZTRM curve with the relaxation curves 
for longer waiting times, using the //-scaling approach 
and taking tff as the scaling parameter. No value of 
\x can give even approximate scaling, especially at short 
times. It turns out, however, that the ZTRM curve can 
be merged with the other TRM curves (using the same 
\x as used for the longer waiting times) , if a much shorter 
characteristic time for the ZTRM is introduced. We call 
this time the "effective scaling time", tff. It is several 
times shorter than the effective cooling time. The inset 
of Fig. 2 exhibits logarithms of and t e ^ for differ- 
ent temperatures. The error bars for the logarithms are 
typically 0.1, but about 0.2 for the lowest temperature. 

The fact that one can consider at least two character- 
istic times for the ZTRM agrees with what is expected 
from relaxation in the positive-temperature-cycle experi- 
ment. The spin-glass behavior at short observation times 
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FIG. 2: The logarithmic relaxation rate for the ZTRM. 
The measurement temperatures, from top to bottom, are 
T/T g = 0.87, 0.77, 0.68, 0.58, 0.5, 0.4, and 0.3. The inset shows 
temperature dependence of logarithms of the effective cooling 
time, tff (solid symbols), and the effective scaling time, tf^ 
(open symbols). 



is governed by t e P , which can be associated with ordinary 
aging effects during the slow approach of the measure- 
ment temperature. However, the long-time relaxation is 
strongly influenced by the metastable states, separated 
by high barriers, resulting from the temperature change. 
Thus, the relaxation rate peaks at tff, and not at tj* . 
This argument can be generalized to include the TRM ex- 
periments with finite waiting times. The logarithmic re- 
laxation rate has a peak at t|y , the effective waiting time. 
It is.. well known that t$ for the TRM is greater than 
t w E3 We believe that this shift in the maximum of the 
relaxation rate is caused by those long-time metastable 
states, which are created during the cooling process. 

It is important to note that memory of the initial state 
depends not only on the cooling protocol, but also on 
the overall complexity of the free-energy landscape. This 
can be illustrated by measurements of the ZTRM in the 
presence of a high constant field. The minimum possi- 
ble overlap, q m in{H), between- two states increases with 
increasing magnetic field, HS3 Therefore, certain con- 
straints on barrier heights are imposed by the field, lead- 
ing to a faster relaxation. The experiment is performed 
in the following way. First, the sample is cooled down 
to the measurement temperature at the field H + AH. 
The field is changed to H, and the zero waiting time 
magnetization (ZTRM) is measured for some time. Sec- 
ond, the sample is warmed up to above T g , and cooled 
down again, all in the presence of the same field H . The 
field-cooled magnetization (MFC) is measured for the 
same time. The difference of these two magnetizations 
is then fitted to a power law at long observation times: 
ZTRM(t, H, AH)- MFC (i, H) cx t~ x W. The values of 
H were varied from 50 Oe to 1000 Oe. The field change 
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FIG. 3: Time dependences of the ZTRM and MFC, measured 
at T/Tg = 0.79 and H = 500 Oe. The field change is AH = 
10 Oe. The inset shows the relaxation exponent, X(H), as a 
function of H 2/3 . 
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FIG. 4: Experimental time dependence of the \trm , tw ) 
and xiRM(t,t w ) for T/T s = 0.87 and t w = 1000 s. The TRM 
relaxation for zero waiting time is referred to as ZTRM. Note 
its large magnitude in comparison with the TRM for t w 7^ 0. 



was small, AH = 10 Oe, so that the response was always 
linear. 

Fig. 3 exhibits the measured ZTRM and MFC decays 
for H = 500 Oe. One can see that time dependence of the 
MFC cannot be neglected at high fields. The inset shows 
the relaxation exponent X(H). The exponent appears to 
be a linear function of q m in oc H 2 / 3 . This is not surpris- 
ing, because its temperature dependence is also close to 
linear, and the AT critical line, T g — T C (H^ oc H 2 ! 3 , has a 
profound effect on spin-glass dynamics EB The estimated 
value of X(H) near the AT line, which corresponds to 
Hat ~ 1400 Oe at this temperature, is 0.27 ± 0.03. It 
is interesting that this number compares well with the 
Monte Carlo result—A w 0.25, for the 3D EA model at 
zero magnetic field. 13 

The observed increase in X(H) suggests that energy 
barriers, created by the cooling process, are lower, if cool- 
ing is performed in the presence of a high field. This 
happens because the accessible phase space is limited by 
Qmin(H). The low-field initial state in our experiments 
appears to be more complex, i.e. characterized by higher 
barriers, than the initial state in Monte Carlo simula- 
tions. 



B. Subaging or superaging? 

It was shown in Sec. II. A that the TRM relaxation can 
be decomposed into a sum of two decays: the IRM, which 
is a response of the system associated with the waiting 
time, t w , and the ZTRM, which is a response of the initial 
state. 

Fig. 4 displays the experimental TRM and ZTRM 
curves, and their difference, for T/T g = 0.87 and t w = 
1000 s. One can see that the response of the initial state 



dominates the measured TRM relaxation. This effect is 
particularly pronounced at low temperatures, where ther- 
mally activated processes are slow. Fig. 4 suggests that 
the ZTRM cannot be neglected, nor can it be regarded 
as a constant. Therefore, the TRM and IRM relaxation 
curves will have different scaling properties, and should 
be studied simultaneously. 

In this paper, we consider a multiplicative representa- 
tion for the measured relaxation curves, in analogy with 
Eqs. (7) and (8): 

XTRM — 4>(t/t w ) PTRA.l(t,t w ,t c ) ; (11) 
XIRM = <f>(t/t w ) pIRM{t,t w ,t c ) ■ (12) 

Here we made an assumption that both the TRM and the 
IRM are characterized by the same scaling part, <fi(t/t w ). 
This is certainly correct in the limit of very long wait- 
ing times, when the two functions are essentially equal. 
Analysis of experimental data, based on Eqs. (11) and 
(12), suggests that this approach works well on the real 
time scales. 

Systematic departures from full t/t w scaling are de- 
scribed by the factors ptrm and pirm- If violations of 
t/t w scaling were only due to the quasiequilibrium behav- 
ior, these factors would depend only on the observation 
time, t. We assume, however, that there is an interplay 
between the cooling and waiting-time effects. Thus, we 
include the waiting time, t w , and an additional argument 
i c , which indicates dependence on the cooling process. 

The correlation function, C(t + t w ,t w ), defined by 
Eq. (1), is independent of t w a t t . — . 0. Numerical studies 
of the SK model demonstrateE§E3 that the correlation 
curves for different waiting times, plotted versus t/t w , 
cross at one point, corresponding to t ~ t w . This is not 
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FIG. 5: Relaxation curves XiRKi{t,t w ) and XTRM{t,t w ) for 
two waiting times, normalized by one at t — t w . The temper- 
ature is T/Tg = 0.87. The IRM decay as a function of t/t w 
is slower for longer waiting times ("superaging"). The TRM 
decay is faster for longer waiting times ( "subaging" ) . 



FIG. 6: Factorization of the relaxation curves from Fig. 5. 
/ = {xirm Xtrm) 1 ^ 2 is an average over two different ini- 
tial conditions, which shows systematic improvement in t/t w 
scaling, g = (xirm /xtrm) 1 ^ 2 describes departures from the 
average and does not scale well. 



surprising, because transition from the quasiequilibrium 
to the aging regime occurs at t « t w , and the correlation 
at this point is near qEA- If the linear response suscep- 
tibility depends on its time arguments only through the 
correlation function, i.e. x = x(C), the susceptibility 
curves should also cross at one point. Experimental re- 
sults do not exhibit this property. The TRM curves for 
long waiting times lie below the curves for short waiting 
times, when plotted vs. t/t w . The IRM curves demon- 
strate the opposite behavior. 

In the present paper, we normalize both the TRM and 
the IRM decays by one at t = t w , and treat them as 
"normal" relaxation functions. This approach has sev- 
eral advantages. First, all departures from full t/t w scal- 
ing, reflected in the shapes of relaxation curves, can be 
observed clearly. Second, the shapes of TRM and IRM 
decays can be compared in detail, regardless of the fact 
that their magnitudes are different. Third, experimen- 
tal results can be directly compared with results of nu- 
merical simulations. We believe that this approach is 
physically justified because the influence of the cooling 
process leads to systematic changes in the shapes of re- 
laxation curves, as discussed in Sec. III. A. This method 
is differ|en.t from the one, traditionally used to study t/t w 
scalingfltl It provides new insights into aging behavior of 
real spin glasses. 

Fig. 5 exhibits the TRM and IRM relaxation curves, 
measured at T/T g = 0.87, for two waiting times, t w — 
1000 s and t w = 6310 s. The TRM data demonstrate the 
familiar subaging pattern: the relaxation, plotted versus 
t/t w , is faster for the longer waiting time. The IRM 
results show a quite different behavior: the relaxation as 
a function of t/t w slows, down as t w increases. This effect 
is called "superaging" 



The phenomenon of superaging may seem rather un- 
usual from an experimentalist's point of view, because it 
has never been observed in TRM experiments. It turns 
out, however, that this is a predominant effect ia Monte 
Carlo simulations. pFig- 4 in Marinari et al.,£3 Fig. ,-4 
in Takayama et aL,ll3 and Fig. 4 in Cugliandolo et a/.EJ 
clearly demonstrate the superaging behavior. In all these 
simulations of the SK model, the correlation function, 
C(t+t w ,t w ), plotted versus t/t w , decays slower for longer 
waiting times. In the case of the EA model, it is more 
difficult to distinguish between the subaging and super- 
aging effects, because the correlation curves exhibit good 
t/t w scaling in the aging regime. However, the relax- 
ation exponent A(T, t w ), defined as C(t + t w ,t w ) oc t~ x 
for t 3> t w , decreases slightbz_with increasing t w , accord- 
ing to Fig. 2 in Kisker et al&3 This fact suggests that the 
superaging behavior appears in the EA model as well. 

Because the experimental TRM and IRM curves have 
the opposite scaling properties (subaging versus super- 
aging), the functions ptrm and pirm, describing depar- 
tures from full t/t w scaling in Eqs. (11) and (12), may be 
expected to have inverse effects on the function <fi(t/t w ). 
Following this observation, we define two auxiliary func- 
tions, / and g. The function / is a geometric mean of 
Xtrm and xirm , and the function g describes departures 
from this mean: 

f(t,t w ,t c ) = (xirmXtrm) 1 ' 2 ; (13) 

g(t,t w ,t c ) = ( X iR M /XTRM) 1/2 ■ (14) 

With this definition of /, the scaling part, cf)(t/t w ), is 
left unchanged, while deviations from full scaling are av- 
eraged out, at least partially. Physically, the function 
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FIG. 7: The scaling parameter fi as a function of tempera- 
ture for different relaxation functions. Note that fi < 1 for the 
TRM ("subaging"), and fj, > 1 for the IRM ("superaging"). 



FIG. 8: Logarithm of the effective waiting time, tfjf/ , as a 
function of temperature. The solid symbols correspond to 
t w = 6310 s, and the open symbols - to t w — 1000 s. Note 



that tt S > t w for the TRM, and t% } < t w for the IRM 
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/ represents an average over two different initial condi- 
tions. Of course, if both xtrm and xirm are normalized 
by one at t = t w , their geometric mean is close to the 
arithmetic mean. The IRM and TRM relaxation curves 
can be factorized as follows: 

Xirm — f g ; XTRM = f/g ■ (15) 

Fig. 6 exhibits the functions / and g for T/T g = 0.87. 
Both of them demonstrate the superaging behavior, but 
the function / shows much better t/t w scaling. The func- 
tion g is responsible for the faster relaxation in the IRM 
experiments, and the slower relaxation in the TRM ex- 
periments. 

C. Comparison of TRM and IRM 

In order to make the above arguments quantitative, we 
analyze the t/t w scaling of relaxation curves using the \i- 
scaling approach, discussed in Sec. II. B. We optimize [i 
to achieve the best possible scaling over a wide range of 
observation times, with particular attention to the aging 
regime, t > t w . Inclusion of the quasiequilibrium decay 
with the exponent a improves scaling at short times, but 
it does not help at longer times. In the present analysis, 
the quasiequilibrium behavior is not taken into account 
explicitly. This gives values of /i a little further from 
unity, but does not affect any conclusions about scaling 
in the aging regime. 

Fig. 7 exhibits values of the parameter n for the TRM, 
IRM, and their combinations. 

In the case of the TRM, \i lies between 0.8 and 0.9 
for relatively high temperatures, but drops visibly at low 
temperatures, T/T g < 0.4. Similar reaujis have been ob- 
tained for various spin-glass samples.EltJQ The error bars 



in our analysis are typically 0.02. At the lowest tem- 
perature, however, they are about 0.1 for the TRM, and 
about 0.05 for the IRM. 

The values of /i for the IRM are greater than unity. 
They are near 1.1, but tend to increase at low tempera- 
tures. According to Fig. 7, there is a certain symmetry 
with respect to fi = 1 between the values of /x for the 
TRM and the IRM. We have also found that, in the case 
of the IRM, the parameter fi is very sensitive to any im- 
perfections of the scaling analysis. For example, if the 
actual observation time is slightly longer, than the time 
t, used in the analysis, the corresponding value of \x will 
seem closer to unity. This issue is important when a 
superconducting magnet is used, because it takes some 
time to warm up the persistent switch, and then to cool 
it down again. The results for the TRM are less sensitive 
to these effects. 

Fig. 7 also displays /x for the functions / and <?, given 
by Eqs. (13) and (14). The value of \i for the function 
/ is about 1.05, while for the function g it is near 1.20. 
Thus, both the TRM and the IRM relaxations can be pre- 
sented as products of two factors, according to Eq. (15). 
The rapidly changing / exhibits fairly good t/t w scaling. 
The relatively slow g describes departures from this scal- 
ing. This function becomes closer to a constant as t w 
increases, suggesting that both the TRM and the IRM 
curves scale better for longer waiting times. At low tem- 
peratures, however, the differences between / and g tend 
to diminish. 

Fig. 8 exhibits values of logarithm of the effective wait- 
ing time, tffi , for different relaxation functions. The re- 
laxation curves were fitted using a 5-order polynomial fit, 
and then differentiated with respect to log 10 (t). The ef- 
fective waiting time, tffl , corresponds to the maximum 
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FIG. 9: The relaxation exponent, X(T,t w ), as a function of 
temperature. The solid symbols correspond to t w — 6310 s, 
and the open symbols - to t w = 1000 s. The open diamonds 
denote A for the experimental ZTRM. The upper and lower 
straight lines are Monte Carlo results for t w = and t w = 
1000 mcs, respectively (from Kisker et ah, Ref. 13). 



of this derivative. The error bars for log 10 (i^' ? ) are typ- 
ically 0.05, but about 0.1 for the lowest temperature. It 
can be seen from Fig. 8, that t^/ > t w for the TRM, 
but t$ < t w for the IRM. The values of the effective 
waiting time for the function / are greater than t w , but 
they all differ from t w by less than 3%. Therefore, the 
average of the TRM and the IRM gives better results for 
the effective waiting time as well. 

Fig. 9 provides a further insight into the nature of 
these phenomena. It displays values of the relaxation 
exponent, X(T 7 t w ), defined through x(Mt«) °t t~ x ( T ' tm ^ 
for t 3> t w . This definition is based on the fact that 
time dependence of spin-glasg-telaxation is essentially al- 
gebraic in the aging regime .LrQlLij The value of X(T,t w ) 
is determined from the slope of the linear fit to log 10 (x) 
as a function of log 10 (t) . Because of time limitations in 
our experiments, the fitting was performed in the inter- 
val t/t w — 1.5. ..10 for t w = 1000 s, and in the interval 
t/t w = 1.5. ..5 for t w = 6310 s. In the case of the ZTRM, 
the effective cooling time, , was used instead of t w . 
The algebraic approximation works well, and the error 
bars for the relaxation exponents are about 1%. 

The first conclusion one can draw from Fig. 9 is that 
the zero waiting time decay, ZTRM, is the slowest re- 
laxation measured in the aging regime. This result is in 
sharp contradiction with results of numerical simulations, 
which show that the zero waiting time decay is the fastest 
possible relaxation. The dotted (upper) line in Fig. 9 rep- 
resents A(T) for t w — from the, Monte Carlo studies of 
the EA model by Kisker et ali3 Our experimental val- 
ues of the A(T) for the ZTRM are about 3 times lower. 
We have also tried to estimate the relaxation exponent 
for the SK model, using numerical results of Marinari et 



aZ.EH and Takayama et a/.El, and fitting them in the in- 
terval t/t w = 1...100. In the first case, for T/T g = 0.5 
and t w = 8 mcs, we find A(T, t w ) ~ 0.15. In the second 
case, for T/T g = 0.4 and t w — 16 mcs, we determine 
X(T,t w ) w 0.20. These numbers are similar to those for 
the EA model. 

Fig. 9 demonstrates that the TRM decay at t > t w 
becomes faster as the waiting time increases, while the 
IRM relaxation becomes slower. The corresponding val- 
ues of A(T, t w ) move towards each other and towards the 
dashed (lower) line, which represents-jiumerical results 
for the EA model with t w = 1000 tocs.E3 One can see that 
the function /, the average of the TRM and the IRM, 
has relaxation exponents, which (at relatively high tem- 
peratures) are very close to those from the Monte Carlo 
simulations. These exponents decrease slightly with in- 
creasing t w , as the numerical results do. 



D. A possible origin of fi 

We are now in a position to offer a tentative explana- 
tion for these phenomena. It is based on two observa- 
tions. First, there is a certain degree of symmetry be- 
tween the TRM and the IRM with respect to the "ideal" 
case, according to Figs. 7-9. Of course, this symmetry is 
not exact and depends on a particular realization of the 
cooling process. Second, when the waiting time, t w , is 
very large, the TRM and IRM become essentially equal, 
according to Eq. (6). This_.situation corresponds to the 
weak long-term memory. Ertll These two observations sug- 
gest that the experimental TRM and IRM decays ap- 
proximate the true spin-glass relaxation from the oppo- 
site sides. 

Let us first consider a typical TRM experiment. When 
the cooling is over, a number of metastable states, sep- 
arated by free-energy barriers of different heights, are 
already populated. We shall refer to them loosely as the 
initial state. The response of this state, the ZTRM, is 
slow and cannot be characterized by a single well-defined 
time scale. The initial state is not random: it has already 
evolved towards some equilibrium state, and this makes it 
energetically favorable. During the waiting time, evolu- 
tion continues in the same direction. Thus, in the case of 
the TRM, aging and cooling effects work together. The 
spin-glass state appears to be older. The characteris- 
tic barrier is higher than A(T, t w ) = ksT\u(± w /tq), the 
barrier, associated with the waiting time t w E3 The TRM 
decay is slower than the ideal t^-dependent relaxation, 
which can be characterized by some exponent Ao . As the 
waiting time increases, influence of the initial state di- 
minishes. The relaxation, plotted versus t/t w , becomes 
faster, and its dependence on t w - stronger. This leads 
to the subaging behavior. Therefore, for the TRM exper- 
iments, we would expect // < 1, t%f* > t w , and A < Ao. 

Let us now discuss a typical IRM experiment. In this 
case, a small magnetic field is applied after the cooling 
process and before the waiting time. The initial state, 
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created by the cooling, is again partially equilibrated. 
The subsequent field change makes it energetically unfa- 
vorable. This is a consequence of chaotic nature of the 
spin-glass state with respect to magnetic field. Even if 
the field change is very small, the new equilibrium state 
(reached at the end of the-,relaxation process) is very 
different from the old one£3 The system can no longer 
evolve in the same direction. Moreover, it has to escape 
from the initial state. Therefore, in the case of the IRM, 
aging and cooling effects work against each other. Dur- 
ing the same waiting time, t w , the system will have to 
overcome barriers, produced by the cooling process, and 
explore a different part of the phase space, thus exhibit- 
ing aging. As a result, the spin-glass state looks younger. 
The characteristic barrier, corresponding to the waiting 
time, t w , is now lower than A(T,t w ) for the ideal aging. 
The IRM decay is faster than the ideal relaxation. How- 
ever, memory of the initial state weakens as the waiting 
time increases. The IRM relaxation as a function of t/t w 
becomes slower and more dependent on t w . This leads 
to the superaging behavior. Therefore, the IRM experi- 
ments should yield /x > 1, t e JJ < t w , and A > Ao- 

If this explanation is reasonable, the observed differ- 
ence in scaling properties of the TRM and the IRM is 
a result of the different initial conditions for aging. De- 
partures from full t/t w scaling are present in the aging 
regime not because something goes, wrong at long waiting 
times (e.g. the interrupted agingEl), but because results 
for short waiting times are strongly affected by the initial 
state. An average over different initial conditions should 
improve t/t w scaling. 

Presence of many metastable states, separated by high 
barriers, is a well-known problem in Monte Carlo simu- 
lations, involving thcrmolization of large samples. The 
system easily gets trapped in a metastable state and 
cannot efficientlyj-eKplore the entire phase space to find 
the ground state.Eil The cooling process in spin-glass ex- 
periments gives a similar result: the system is trapped 
before the waiting time begins. However, Monte Carlo 
studies of aging phenomena do not simulate the cool- 
ing protocol. The simulations are started directly at the 
measurement temperature from a random initial config- 
uration, a nd , r es ult s are averaged over different initial 
conditions .E3o'Ej'E3EZI This may be the reason why the 
subaging behavior is not normally observed in numerical 
simulations. Moreover, the random initial configuration 
corresponds to zero net magnetization. This means that 
Monte Carlo studies are probably closer to the IRM, than 
to the TRM, experiments.liZi The difference between the 
two is not very important, if memory of the initial state 
disappears quickly. This, however, is not the case in real 
experiments. 

The factorization of the TRM and IRM relaxation 
curves, Eq. (15), allows us to illustrate a possible origin 
of /i in a simple way. The functions / and g, Eqs. (13) 
and (14), can be approximated algebraically in the aging 
regime: 

/ K (i/er Ai ; 5 «(t/er A2 . (ie) 
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FIG. 10: Illustration of a possible origin of /i. Parameters 
/xi and Ai correspond to / = (xirm Xtrm) 1 ^ 2 , while /12 and 
A2 correspond to g = (xirm /xtrm) 1/2 ■ 

It follows from Eq. (15) that the relaxation exponent for 
the IRM is Ai + A2, while the exponent for the TRM is 
Ai — A2. The corresponding values of /i depend on the 
ratio r = A2/A1: 

M/.RM = (Ml + ^2 r)/ (1 + r) ; (17) 

Vtrm = (Mi - r)j (1 - r) . (18) 

The experimental values of \i\ and \ii are exhibited in 
Fig. 7. The exponents Ai are displayed in Fig. 9, and the 
ratios r = A2/A1 - in the inset of Fig. 10. 

The values of n for the IRM and the TRM, determined 
from Eqs. (17) and (18), are shown in the main body of 
Fig. 10. They compare well with the experimental val- 
ues of /x, obtained from analysis of the actual relaxation 
curves. This comparison demonstrates how the func- 
tion g, which exhibits strong superaging properties with 
(i 2 ~ 1.2, leads to fi ~ 0.9 for the TRM, and fi ~ 1.1 
for the IRM. The results in Fig. 10 were obtained us- 
ing values of r for t w = 1000 s. The values of r for 
t w = 6310 s give similar numbers. For the sake of sim- 
plicity, the observation time, t, is used in Eq. (16) instead 
of the /u-dependent effective time, t e , given by Eq. (9). 
The fact that Eqs. (17) and (18) yield values of \i close 
to the their actual values suggests that this approach is 
justified. 

Appearance of the /x-scaling can be explained in the 
following way. Imagine that there is a function that ex- 
hibits good t/t w scaling. Let us also assume that effect of 
initial conditions is described by another function, which 
does not scale well. This second function will necessarily 
display superaging behavior, approaching a constant as 
t w increases. Imagine also that both these functions can 
be approximated by power laws of time, which is often 
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the case in spin-glass dynamics. If the first function is 
fast, and the second function is slow, the resulting decay 
will exhibit the /z-scaling with fi not far from unity. Of 
course, this is only one of possible scenarios. 

According to Eqs. (17) and (18), in order to achieve 
the best scaling of both the TRM and the IRM (with 
fi close to unity), one has to make the ratio r = X2/X1 
as small as possible. The inset of Fig. 10 suggests that 
r decreases with temperature and can be well approxi- 
mated by a linear function of T/T g . It also diminishes 
with increasing t w . Thus, higher temperatures and longer 
waiting times should improve t/t w scaling of measured 
relaxation curves. 

Eq. (18) also predicts that the quality of /i-scaling of 
the TRM curves should become worse as r — > 1 and 
fii — > fX2, because of the uncertainty 0/0. This happens 
at low temperatures: the experimental values of /j, for the 
TRM drop sharply and their error bars increase consid- 
erably. The physical reason for this is simple. The TRM 
relaxation at low temperatures is completely dominated 
by the ZTRM, which is very slow. The waiting time ef- 
fects become almost indistinguishable, and the t/t w scal- 
ing does not make much sense anymore. In this case, the 
TRM is nearly constant, and both / and g, Eqs. (13) and 
(14), are governed by the IRM. Therefore, the simple av- 
erage of the TRM and IRM does not give much better 
scaling at low temperatures, as seen in Fig. 7. 

These results can be interpreted in terms of an effec- 
tive change in barrier heights. Let us consider the loga- 
rithmic relaxation rates: Str.m = — dXTRM /d\~n.(t) and 
Sirm = —dxiRM I '91n(f). As usually, we assume that 
both xtrm and xirm are normalized by one at t = t w 
before differentiation is performed. A ratio of Strm and 
Sirm can be regarded as a probability of overcoming a 
barrier: 



(Ixtrm/ dxiRM = exp(-A /k B T) 



(19) 



This formula defines the effective barrier Aq(£, t w ), de- 
scribing the difference in relaxation rates for two different 
initial conditions. 

Fig. 11 exhibits values of Ao/fc^T at t — t w , extracted 
from the experimental data according to Eq. (19). One 
can see that they grow steadily as temperature decreases. 
They are also comparable with differences between the 
values of ln(t^) for the TRM and the IRM. Values of 
Ao(t w ,t w ) are displayed in the inset of Fig. 11. They 
change only slightly over a wide range of temperatures, 
and decrease with increasing t w . The error bars for Ao 
are about 0.05. 

It has been shown experimentally that a change in a 
barrier height, 8 A, under a .temperature change, 8T, is 
independent of temperature E3 The same appears to be 
true for Ao, which corresponds to an average change in 
barrier heights due to the change in initial conditions. 
This fact suggests that Ao is a result of cooling. Both the 
TRM and the IRM are influenced by the cooling process, 
and it is mainly this influence, that leads to systematic 
deviations from full t/t w scaling in the aging regime. As 
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FIG. 11: The effective barrier Ao, describing difference in 
relaxation rates for the TRM and the IRM, at t = t w . Note 
that it is almost independent of temperature, as shown in the 
inset. 



the waiting time increases, Ao diminishes, which means 
that memory of the initial state is gradually washed away 
by the aging process. 

It has been already mentioned that evolution of the 
system during the waiting time t w can be characterized 
by the effective barrier A(T,t w ) = fc^Tln^/Vo). The 
barriers are raised by Atrm for the TRM, and reduced 
by Airm for the IRM. The average difference in barrier 
heights for these experimental protocols is A = Atrm + 
Airm- Using Eq. (19), one can determine A only. Fig. 8 
suggests that Atrm and Airm may be different, but 
have similar properties. Therefore, what is true for Ao, 
is likely to be true for Atrm and Airm taken separately. 

The characteristic barrier, A(T,t w ), decreases with 
decreasing temperature, because it corresponds to 
thermally activated aging dynamics. The barrier Ao is 
almost independent of temperature, presumably because 
it is related to the cooling process. Therefore, departures 
from full t/t w scaling should become more pronounced 
as temperature is lowered. Even if the cooling is fast, 
yielding small Ao, it will always be possible to find a 
sufficiently low measurement temperature, T, so that 
A(T,t w ) ~ Ao, and the t/t w scaling is strongly violated. 
According to this argument, the sharp drop in fi for the 
TRM at low temperatures is a natural consequence of 
cooling effects. 

Some of the measurements, discussed above, have been 
repeated at CEA Saclay. A SQUID susceptometer S600 
by Cryogenics Ltd (UK) was used for this purpose. The 
cooling process, which is a unique characteristic of an ap- 
paratus, was somewhat different from the one exhibited 
in Fig. 1 . The experimental data, however, are consistent 
with the results, reported in this paper. 
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IV. CONCLUSION 

There are many different approaches to analysis of t/t w 
scaling in spin glasses. In this paper, one more method 
is proposed. It is based on the idea that the departures 
from t/t w scaling, observed experimentally in the aging 
regime, may reflect the influence of the cooling process. 
This means that the thermal history of the spin-glass 
state cannot be neglected in analysis of aging phenomena. 
The initial condition for aging is not random in spin- 
glass experiments, and has a profound effect on scaling 
properties of measured relaxation. It would be incorrect 
to say that the cooling effects are spurious. They are as 
important, as the aging phenomena. This is because thejz. 
reflect the hierarchical nature of spin-glass dynamics,E£| 
revealed by temperature-variation experiments. 

The approach, used in this paper, is based on the de- 
tailed comparison of time scaling properties of the TRM 
and the IRM. We argue that, because of strong memory 
of the initial state, neither of these functions represents 
the true i^-dependent relaxation. Even though they 
are characterized by the same waiting time, their t/t w 
scaling properties are remarkably different. We exploit 
this difference by taking an average of Xtrm and Xirm- 
This averaging over different initial conditions leads to 
systematic improvement in t/t w scaling. Factorization 
of measured relaxation functions allows us to explain 
scaling properties of both the TRM and the IRM in 
a unique way. Our results suggest that the observed 
deviations from t/t w scaling in the aging regime ara 
largely due to cooling effects. The interrupted aging,B 
which is another possible reason for the lack of scaling, 
may be expected to play a role at longer waiting times, 
especially in polycrystalline samples. 



The conclusions, reached in this paper, are by no 
means final. Their verification requires experiments 
in which different cooling protocols in the immediate 
vicinity of the measurement temperature can be imple- 
mented. Another option is a rapid field quench from 
above the AT line. Monte Carlo studies, including 
simulations of the cooling process in addition to aging 
effects, may also be very helpful. The problem of 
t/t w scaling in spin-glass dynamics deserves further 
experimental and theoretical investigation. 

NOTE. After this paper had been completed, we re- 
ceived a preprint by Berthier and Bouchaud,E3 devoted 
to numerical studies of the Edwards- Anderson model in 
3 and 4 dimensions. One of their conclusions is that "a 
finite cooling rate effect... leads to an apparent sub-aging 
behavior for the correlation function, instead of the 
super-aging that holds for an infinitely fast quench" . 
This is in agreement with the conclusion of the present 
paper that t/t w scaling properties strongly depend on 
the initial condition. 
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